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Abstract. Radiative corrections to the single top s and t channel production processes are revis- 
ited. Complete one-loop electroweak corrections are calculated within the SANC system. New is a 
study of the regularisation of the top-legs associated infrared divergences with aid of the complex 
mass of the top quark. A comparison of these electroweak corrections with those computed by 
the conventional method is presented both for top production and decays. Standard FORM and 
FORTRAN SANC modules are created. These modules are compiled into a package sanc_cc_vl . 40, 
which may be downloaded from SANC project homepages. Numerous numerical results are pre- 
sented at the partonic level with the aim to demonstrate the correct working of modules. These 
modules are intended to be used in Monte Carlo generators for single top production at the LHC. 
Where possible, we compare our results with those existing in the literature; in particular, a 
comprehensive comparison with results of the CompHEP system is given. 
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1 Introduction generators, such as MC@NLO [TU] or SingleTop [TT] . 

incorporated in the standard LHC tools. 



In continuation of work on electroweak one-loop cor- 
rections to processes involving the top quark [TJ[1|, 
we present here a calculation of EW one-loop correc- 
tions to single top quark production and top quark de- 
cay. To incorporate these processes in the SANC frame- 
work [3] is a natural extension of the previous work. 

The interest in the study of single top production 
has been stated many times (see e.g. [1] and references 
therein), suffice it to say that this is the only way of 
measuring the CKM matrix clement \Vtb\ directly and 
thereby providing a sensitive test of the 3-generation 
scheme of the Standard Model. Indeed, a significant 
deviation of \Vtb\ from the value demanded by uni- 

tarity of the CKM matrix would be an indication of In tllis P a P er we calculate the electroweak correc- 
the presence of a fourth generation of fundamental tions to the 3 and t channels of single top quark pro- 

c duction and revisit the top quark decay El- The im- 

termions. 1— ' 



There are three channels of single top quark pro- 
duction: t channel, s channel and associated produc- 
tion. The t channel process, bq — > tq' , where q and 
q' are light quarks, is the channel known to have the 
largest cross section [T2HT5] . Calculations for the LHC 
at 14 TeV have led to the prediction of this cross sec- 
tion being more than an order of magnitude greater 
than that of the s channel process qq' — > tb and about 
four times greater than the cross section of associated 
production qg — > tW (see e.g. [IS]), and it has been 
stated that the electroweak corrections are of compa- 
rable magnitude to the QCD corrections [12] . 



portant new feature of the present work is the inclu- 
sion in the calculations of the width r t of the top 
quark. This is a nontrivial step, and we will show that 
the results are sensitive to r t . 



Precision calculations of single top quark produc- 
tion have been done for a number of years (see e.g. [5] 
and references therein) but remain of continued inter- 
est. This is motivated on the one hand by the obser- 
vation of such events at the Tcvatron DO 0(7] and The paper is organized as follows. 
CDF [8K] and, on the other hand, by the need to pre- l n section 2 we present the covariant amplitudes for 
pare software to analyse single top quark events at the the processes under consideration, working within the 
LHC, running at 7 TeV. framework of the multichannel approach (see R,ef |17j ). 

Most of the theoretical work in single top produc- In section 3 we introduce the infrared rcgularisation 
tion has been concerned with higher order QCD cor- by the complex mass of the top quark, at first for the 
rcctions, leading to the development of Monte Carlo virtual QED corrections, and then for the hard photon 
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radiation. Numerical results are presented in section 
4 and the Conclusions and Outlook in section 5. 

2 Covariant Amplitude 

In this section we proceed in the spirit of "multi- 
channel approach", see Ref. [17] ■ First, we consider 
annihilation into vacuum with all particles incoming. 

2.1 All particles incoming 

It is convenient to consider two vacuum diagrams: 



Covariant Amplitudes (CA) of Figs. Q] and [2] arc 
characterized by four different structures and scalar 
Form Factors (FF) if the mass of the light quarks is 
neglected the and b quark mass is not neglected. Omit- 
ting Dirac spinors, one may write a common expres- 
sion for this CA of the process in terms of scalar form 
factors, J 7 LL (s,t), see [3]: 




t w 

Fig. 1. The tbud — > processes 



and 




A = i e 



2 d w (s) 



7m (! + 7s) <8> 7m (1 + 75) Fll{s, t) 
+7m (1 - 75) 8> 7m (1 + 75) ?hl{s, t) 
+ (1 + 75) 8> 7m (1 + 75) {-iD^nis, t) 
+ (1 - 75) ® 7m (1 + 75) {-iD tl )F RD {s, t)] , (1) 



with 



-Dm = (Pi - Pa) 



(2) 



and 4-momentum conservation reads 



Pi + P2 + P3 + Pi = , 



the invariants are defined by 



(3) 



-(Pl+P2) 2 , t = -(p 2 +p 3 )' 2 



and 



1 



(4) 



(5) 



Fig. 2. The tbud — > processes 

Here black ovals symbolically denote all one-loop 
insertions to the corresponding tree diagrams. 



2s 2 w s - Ml + iM w r w ' 
Here V = VtbV u d is the relevant product of CKM ma- 
trix elements and the symbol ® means, considering 
the factor of Tll in diagram Fig. [TJ 

7m 7+ ® 7m 7+ = u(p 2 )^ tl j+v(p 1 )v(p 3 )-f^ + u(p 4 ) 
with 7+ = 1 + 75 
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and appropriate changes in the other cases. The scalar 
form factors T are labeled according to their struc- 
tures, see [5]. 

2.2 Conversion to top decay 

The CA for the decay t(p 2 ) -> b(pi) + u(p 3 ) + d(pi) 
is derived from Fig. [T] by the following 4-momentum 
replacement: 

Pi -Pi, P3 -> -Pa, 
P2 -> P2, Pi -Pi- 
For the decay t{p\) — > b(p 2 ) + d(p 3 ) + u(jp±) it is 
more transparent to make the replacement from Fig. [2] 

Pl -> pi, p 3 ^ ~P3, 

P2 -> ~P2, Pi -> ~Pi- 

As a result one has two decay diagrams, schemat- 
ically representing one-loop EW corrections: 




yy nn. ps 

t u 

Fig. 3. The t — > bud processes 

From the diagram of Fig.|3]thc four Helicity Ampli- 
tudes (HA) can be derived by the standard techniques 
used in SANC, sec Eqs. (47-50) of Ref. 0. 

For the case of the process of diagram of Fig.[3]the 
HA are similar, although not identical. For their exact 
expressions sec the relevant module in the SANC tree. 




2.3 Conversion to s channel 

The CA for the s channel single-top production pro- 
cesses u(pi) + d(p 2 ) — > b(ps) + i(pi) is obtained from 
Fig. [U by the following 4-momentum replacement: 

Pl -> ~P3, P3 -> Pl, 

P2 -> -Pi, Pi -> P2, 

leading to the diagram 




Fig. 5. The ud — > tb process 

For the processes d{pi) + u(p 2 ) — > t(p 3 ) + b(pi) from 
Fig. [5] the conversion 

Pl -> ~P3, P3 -> Pl, 

P2 -> ~Pi, Pi -> Pi, 

leads to the diagram of Fig. [51 

For the corresponding HA one can get rather com- 
pact expressions. 

• Helicity amplitude for ud — > bi 

%+ = k sm${(m t ki + m h k 2 )T LL 
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Fig. 6. The du — > tb processes 

+(m t k 2 + m b ki)F L R 
-P + P-(k 2 F LD -kiFR D )}, 

1L+-++ = k sind{{m t k 2 + mbki)J r LL 
+(m t ki + m b k 2 )F L R 
+P + p-{k 1 F LD -k 2 F RD )}, 

H+ — + = ko\fsc + (kxF L L + k 2 F L n), 
= koy/sc— (k 2 J-LL + kiF LR ). 

Hclicity amplitude for du — > tb 

%+ = k sinz9{(TO t fc 2 + m b kx)T LL 

+(m t ki + rn b k 2 )F L R 
-P + P-(k 2 F LD -kiF RD )}, 

T-L+-++ = k s'm{){(mtki + m b k 2 )Thh 
+ {m t k 2 + m b ki)F L R 
+P + p-{k l F LD -k 2 F RD )}, 

U+ — i — = fcoVsc_(fc 2 Jii + kiT LR ), 

H+ — + = ko^/sc+^iTLL + k 2 F LR ). 



Here 



ka 



x(M 2 w ,s) 



(G) 



(7) 



h, 2 = P-±P + , 

P = \J s - (m t ± m b ) 2 , 
c± = l±cosi? (8) 

and •& is always the angle Z(pi,p3). 

2.4 Conversion to t channel processes 

In the CAs for the t channel single-top production 
processes b(pi) + u(p 2 ) -> d(p 3 ) + t(pi) and b(pi) + 
d{p 2 ) — > u(p 3 ) + i(j>/i) it is convenient to make the 
replacement "in pairs". From Fig. [1] one may perform 
two 4-momentum replacements: 



Pi 

P2 
P3 



Pi, 
P2, 



Pi 
P2 
P3 
P4 



Pi, 
~P4, 
-P3, 
P2, 



P4 -> -P3, 

which give rise to two different physical t channel pro- 
cesses, described by two symbolic diagrams: 

^ Pi P4 / y , \ Pi Pa 

b /* t b Y \ X 




P2 P3 




P2 P3 



2s 2 w s-M^ + iM wlw 



Fig. 7. The bu — > td and bd — > tu processes 

For the processes b(pi) + u(p 2 ) — > t(p4) + d(p 3 ) and 
b(pi ) + d(p 2 ) — !• t (pA) + u(p3 ) , the pair of replacements 
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from Fig. H 



Pi -> ~Pa, 
P2 -> Pi, 
P3 -> ~P3i 
P4 P2, 



Pi 

P2 
P3 

Pa 



-Pa, 
Pi, 

P2, 
~P3, 



rk 2 

•> RL 

S 

m t k 2 



Frd \, 



fc sintf -N+ 2s tb T LL 



-m t 



gives the corresponding pair of symbolic diagrams for 
the two remaining t channel processes: 



-iV Zmt+bFLL 
m t 



— Sb-F RL - k 2 T L D) - ki^RD | 



Pi Pa 



Pi Pa 




k c+{N+^(^-T RL 

— kiJ-LD t ——J'RD 



Fig. 8. The bu — > td and bd — > tu processes 



-N c_.s b [T RL + rat (T LD - T RD )] } . (9) 

Specific to Eq. ([9]) is the following notation: 

ko = V X (M^,s)N 2 N^ 3V 
ki = 2s t b - c+s b , 

k 2 = 2s tb m b + c + m t Sb. (10) 



The HAs arc all different for the four types of pro- The other notation common to Eqs. © and (flU)) is 
cesses. They are listed below in another sequence than 



they were presented diagrammatically above. First, we 
give the more complicated HAs for 4-particle and 4- 
antiparticle processes: 

• Helicity amplitudes for ub — > td 



s t b = s + m t m b , 
St = s — m t 
c± = 1 ± cosi?, 



mt+b 

Sb 



m t + ra b , 



(11) 



and 



n = fco{iV + yic 

mt 



2m t+ bJ : LL 
— (c-Sb^RL - k 2 T LU ) - kxFRD 

-J\Tc_ \2stbTLL 
h 



N t 



N 

(2,3) 

N+ 



-m t 



/ * 


V s tb c - - 


f sm 2 c+ 


N (3,2) = 


js t s b 
]/ — 


m t+ b 
V2 ' 


N~ = 



Stb 



2s 



(12) 



H 



fco sin^|iV + c + ,S6 



•F RL 



-mt (Fld - Frd) 



H 



• Helicity amplitudes for bu —> td 

fc 2 



= fc |-A+\/ic + 



-Frl 
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(2ro t+6 m t + c_s t ) Tld 



Trd 



-AT~c_ 



2stbJ~ LL 



-mb^—FRL + c + s t (T L D - Trd)^ }, 



H- 



= k sin${N + 



-nib 



-AT 



Zstb^LL 
T RL + C+Sti^LD — Trd) 



rfc 2 



-Trl 



-Kl-TLD -TRD 



}■ 



fc sin # j -iV + c + s t [j" flL 



-AT 



2m t+6 J r LL 

+—(c-S t TRL - k 2 T L D) - kiF RD 

= k a c + {-N + ^ 
m 6 



-St J" RL + Vld) + ki^RD 

-Af~c_s t [J"^ + m b (T L D - Tad)] }(.13) 
Notation specific to Eq. (|T5j) is: 

fci = 2,s tb - c + s t , 

fc 2 = 2s tb m t + c + m b s t . (14) 

Next we give HAs for the i channel process of anni- 
hilation of different flavors; they are simpler than the 
previous ones. 

• Hclicity amplitudes for bd — > tu 

H-++- = fe c + (Ji£ - m t T L D - m b TR D ), 

U-+++ = -fcoVssintf (— Tll 
_ m t m b \ 

J~ LD J~ RD , 



-ko 



m b 
m t m b C- 



T-L+++- = fcoVssini? — (T L l - m t F LD ) - T RD 
L s 



J~LL — 2Trl 



+c + (m b T L D + m t T RD ) 
• Hclicity amplitudes for bd — > tu 



(15) 



U. 



U 



= -k 



mtlTLbC- 



— 2Trl 



+c + (mtTLD + rribFRD) 

-fcov^sini? {—Tll 
\ s 



_ mtmb 

-J~LD -TRD 



H-\ = koy/s siwd 



— (Tll - mbTLD) - Trd 
s 



H+ — + = k c + (T L L - m b T L D - m t T RD ). (16) 
In Eqs. mijl and (HD 

k =V X (Ml,s)P + (s,0,m b )P + (s,0,m t ) , (17) 
with a function typical to hclicity amplitudes: 

P ± (s,x,y) = ^s-(x±yf. (18) 



3 Infrared regularization by the complex 
top quark mass. 

Infrared divergences associated with interactions of 
photons with on-mass-shell top legs in the limit of 
zero width of the top quark may be regularized by 
any conventional method (photon mass, dimensional 
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rcgularisation), but physically they are naturally reg- 
ularized by the presence of the finite top quark width 
r t . Here this approach will be exploited. 

3.1 Top quark width and virtual QED correction 

There are three diagrams contributing to virtual one- 
loop QED corrections involving the top quark: see 

Figs. El HE and QU 




t u 

Fig. 9. The self-energy photonic QED correction for 
tbud — > processes 





t u 

Fig. 11. The direct box photonic QED correction 

and the expression for ^^f D was explicitly expressed 
in terms of standard Passarino-Veltman (PV) func- 
tions Co- To regulate the amplitude the photon was 
given a mass A. However, one can alternatively intro- 
duce the top quark width, i.e. instead of 

C {-ml, -ml, Q 2 ; m t , A, m b ) 



we now write 



C (— m 2 , —ml, Q 2 ; m t , 0, m b ) 



(20) 



with m t = m t + A t , A t = —imtr t . 

The corresponding replacement for the IRD term 
due to Fig. El reads 

-im t r t 



~ 2 

In 

mf 



- 9 
in — 



(21) 



Fig. 10. The vertex photonic QED correction 

The appropriate formalism was developed in Ref. [2] . 
In particular, it was shown that of the four scalar form 
factors JF only Fll contains IRD QED as well as weak 
contributions which can be separated thus: 



H \ "H 

The expression for the amplitude of Fig. [TT] con- 
tains an infrared and mass singular Dq function. In 
Ref. [18] we have shown that these functions can be 
reduced by standard Passarino-Veltman reduction to 
an infrared finite mass singular auxiliary function J aw 
and a mass singular Co function. This has become 
standard SANC approach which we apply here. 

The relevant Co function is 



Fll = 1 



e - QED .9 



— T 
16^ 2 LL 



weak 



16-7T 



(19) Cl R -Co(-m 2 ,-m 2 ,Q 2 ;m t ,0,m 6 ) 



(22) 



s d 









-uJ±) 


V VdlJ 


\VdiJ 



-Li 2 
-Li 2 
-ln( 
-Li 2 
-Li 2 
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3.2 Top quark width and hard photon radiation 

Wc must also consider the amplitude with real photon 
emission off the t quark, see Fig. [T2] This amplitude 
involves a propagator with denominator 



1 - Vdi 
Udi - Vdi 

1 - Vdi 
Vdi - Vdi 



+ Li 2 
+ Li 2 I - 



Vdi 



Vdi - Vdi 

Vdi 
Dd2 - Vdi 







('--) 




V Vd2j 


\VdzJ 



f 1 - Vdi 
\Vd\ - Vd2 

1 - Vd2 
Vd2 - Vd2 



- Li 2 - 



Vd2 



Li 2 



Vdi - Dd2 

Vd2 
Vd2 - Vd2 




where 



B d = Q A + mf - mi + Q A /\Q 2 \ie, 

S d = \J (<3 2 ) 2 + 2Q 2 (m 2 + m 2) + (m 2 _ m 2 )2; 



Vdi 



Bd — Sd 
2Q 2 



Vd2 



Bd + Sd 
2Q 2 ' 



(23) 



and 



Bi=Q 2 + m 2 - m 



S, 
Vdi 



s]{Q 2 ) 2 + 2Q^(m 2 ' ^ - {rh 2 _ m 2 b) 

Bi — Si Bi + Si 



t u 

Fig. 12. The top Bremsstrahlung diagram 

(p2 — P5) 2 + ml — i m t r t which evaluates to z 2 — i m t T t 
with z 2 = -2p 2 ■ P5- 

The amplitude of the hard photon radiation pro- 
cess without taking account of the top quark width 
can be presented in the following form: 



= A^ +A udb w- (27) 



2Q 2 ' 
Its limit for mj — >■ is 



Vd2 



1 *m(±W n 



hi 



Bi 
B d 



2Q 2 



( ^ ) 



(24) 



A, 



Here A; ^ is the amplitude of the radiating top quark 
and Audbw is the rest. The amplitude (|27p is gauge 
invariant: Ay ■ p§ =0. 

One can define the amplitude of the hard pho- 
ton radiation process taking account of the top quark 
width as 



hi 



B, 



In 



B n 



Here 



<t - Li2 (l> 



Bi = Q 2 + m 2 t , 



(25) 



Ay = A t + A u dbw 



-Af> + Audbw,(W) 



B d = Q 



(26) 



z 2 - im t r t 

where the top quark propagator with width _T t was 
introduced. 

The amplitude ([25)1 destroys gauge invariance: 
Ay • P5 7^ 0. The gauge violation is of order r t /mt- 
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To recover gauge invariance one can modify the 
amplitude (|28p in the following way: 



A 



Z2 



-2 



z 2 - im t r t 

(0) 



A 



<lbW 



(29) 



z 2 - im t r t 

Such transformation (called overall scheme, (o)) 
obviously leads to gauge invariance and the mod-squared 
amplitude will have the form 



1 A udbW + A{ 0) (AudtwY 
+ \Aud bW \ 2 - (32) 
In the next chapter we present numerical results 
for all width schemes introduced above in order to 
estimate the numerical impact of gauge violation in 
fixed schemes. 



a over 
I ^7 



4 + ™> 2 t r t 2 



A 



(0) 



+ 



4 0) ) a udbw + A t (AudbwY 



\A U 



<lb\Y 



(30) 



As is seen from equation (|30[) . all terms are now 
regularized by the top quark width but the last term 
is multiplied by an artificial factor which cannot be 
treated in the soft (Born-like) kinematics. 

To avoid the problem with integration one can in- 
troduce a new gauge non-invariant amplitude (which 
we call fixed\ scheme (fif): 

The amplitude A 1 with top quark width intro- 
duced ah initio (|28[) is called fixed\ : 



|^/i:re(ii | 



,2 
-2 

2Re 



Z2 



4 0) 



\Au 



z 2 - im t r t 

2 



A[ (AudbwT 



dbW\ 



(31) 



An alternative version, called fixed 2 , obtains by setting 



equal to one the coefficient of \A U 



dbW\ 



in Eq. ((50]) : 



Afixed 2 



A 



(0) 



4 Numerical Results 

In this section we present the SANC results for the cross 
sections of single top quark production and for the de- 
cay widths of the top quark. The tree level contribu- 
tions, both Born and single real photon emission, are 
compared with CompHEP |19j . All numerical results for 
this section were produced with the standard SANC 
INPUT working in the a(0) EW scheme, if not stated 
otherwise. 

The standard SANC input parameters set reads: 





= 1.16637 x KT 5 GeV~ 2 






a(0) 


= 1/137.035999, 




= 0.00729735257, 


M w 


= 80.403 GeV, 


r w 


= 2.141 GeV, 


M z 


= 91.1876 GeV, 


r z 


= 2.4952 GeV, 


Mh 


= 120 GeV, 


m t 


= 174.2 GeV, 


m u 


= 62MeV, 


m d 


= 83MeV, 


m c 


= 1.5 GeV, 


m s 


= 215 MeV, 


nil 


= 174.2 GeV, 


m b 


= 4.7 GeV, 


\v ud \ 


= 1, 


\Vcs\ 


= 1, 


\v us \ 


= 0, 


\Vcd\ 


= o, 


\v tb \ 


= 1. 




(33) 
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This section is subdivided into three subsections: for 
the decay channel, s channel and t channel; for the lat- 
ter subprocesses results for both direct (bu — > td) and 
crossed (bd — > tu) channels are presented separately. 

Each subsection, in turn, consists of two paragraphs. 
The first paragraphs contain a comparison with Comp- 
HEP of the Born cross section and of the hard photon 
contribution with a cut on photon energy E 1 in the 
rest frame of the top (for decay channel) and in the 
cms of the initial particles for s and t channels. 

In the second paragraphs we present a study of 
complete EWRC at the partonic level. We show the 
decay widths for the decay channel and cross sections 
for the reactions at the tree and one-loop levels. The 
one-loop EWRC is defined as usual: 

^jlloop ^jtrcc 



(34) 



where O is either width (r) or cross section (a). 

Furthermore, the dimensionless soft-hard separa- 
tor Co subdivides soft and hard real photon contribu- 
tions. Their sum does not depend on Co. 

For decays a — > b + c + d + 7 with s = — (p a — Pb) 2 
and s' = M? d and for the processes, a + b^c + d + 'y 
with 



= M£ 6 and s' = Mi d 



Q=l- 



and E~, > —10 . 
7 ~ 2 



(35) 



For each of the four channels we present Tables illus- 
trating the existence of an E y (Co) stability plateau 
for two options r t = 0, 7^ 0, and answering the ques- 
tion how big the difference is of the EWRC between 



the two r t options. For reactions, we also investigate 
the question of the initial quark masses {m q ) indepen- 
dence of the subtracted quantities. 

For the explanation of the necessity of the subtrac- 
tion procedure and its realization in the MS scheme 
within SANC we refer to our earlier papers [ID]- [2Tj . 

4.1 Decay Channel t -> b + u + d 

4.1.1 SANC-CompHEP comparison 

No new comparison for this channel was done com- 
pared to Rcf . 2 . If the top quark width is taken into 
account, SANC and CompHEP numbers agree well. 

4.1.2 One-loop EW corrections 

Process t — > b + /i + + Vn . 



r Born =0. 1490949(2) GeV 




r f 


= 




E 7 , GeV 


HT 1 


10~ 2 


10~ 3 


^-illoop 


0.159503(1) 


0.159495(2) 


0.159499(7) 


S 


6.981(1) 


6.975(1) 


6.978(5) 


-Tf 7^ 0, fixedi 


^lloop 


0.160736(2) 


0.160784(2) 


0.160787(8) 


5 


7.810(1) 


7.841(2) 


7.842(5) 



Table 1. The total lowest-order and one- loop corrected 
widths r Born and T 11 in GeV and relative one-loop correc- 
tion 8 in % for E-, = 10" 1 , 10" 2 , 10" 3 GeV. 
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As is seen from Table [TJ there is good ui stability 
for Erf < 10~ 2 GeV and about +0.9% increase of the 
EWRC for this partial decay width if one uses r t 7^ 
option as compared to the r t — one. 
Process t — > b + u + d . 



r Born =0.4472847(7) GeV 




r, 


= 




E^, GeV 


10- 1 


10~ 2 


io- :i 


^-illoop 


0.47922(1) 


0.47920(1) 


0.47918(1) 


S 


7.139(1) 


7.135(1) 


7.131(3) 


r t 7^ 0, fixedi 


^-nlloop 


0.48293(1) 


0.48312(1) 


0.48311(1) 


s 


7.969(1) 


8.012(1) 


8.010(2) 



VS/GeV 


ws 


200 


1000 


7000 


r t = 


CompHEP 




8.172(1) 


17.19(1) 


0.7966(2) 


SANC 




8.173(1) 


17.19(1) 


0.7964(1) 


r t = 1.54688 GeV 


CompHEP 





7.835(1) 


16.93(1) 


0.7882(1) 




f 


7.788(1) 


16.90(1) 


0.7884(1) 


SANC 





7.834(1) 


16.93(1) 


0.7885(1) 




fl 


7.788(1) 


16.91(1) 


0.7884(1) 






7.787(1) 


16.91(1) 


0.7884(1) 



Table 2. The same but for t — > b + u + d process. 

As is seen from Table [2j one can draw similar con- 
clusions for both considered partial decay channels. 

4.2 s channel 

There are two generic s channel processes: 

d + u^b + t, and u + d -> t + b. (36) 

4.2.1 SANC-CompHEP comparison 

The numbers of this subsection are produced with 
It is sufficient to consider only one: u + d —> t + b. For SANC setup, Eq. (|33[) . The aim is to demonstrate the 
this comparison we use the CompHEP (v. 4. 5. 1) setup intervals of stability of one- loop corrected EW cross 
with a cut on the cms photon energy _E 7 > 2 GeV. sections cr Uoop and relative EWRC S against variation 



Table 3. Comparison of the cross section a ha 
(fb) for three cms energies, two options: Ft = 0, 7^ 0, and 
three width schemes (ws): o-overall, f-fixed (ti(2) = fi xe di(2) 
for the case of SANC). 

As is seen from Table [31 there is good agreement 
of numbers obtained from SANC and CompHEP within 
the statistical errors for all considered options and cms 
energies. It is worth to note that the two versions of 
width schemes agree well for all energies except near 
the threshold and that the fixedi scheme agrees with 
the CompHEP fixed scheme as expected. 

4.2.2 One-loop EW corrections 
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of the soft-hard separation parameter, Q and to study 
the difference between the two options: r t = 0, 7^ 0. 
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5=200 GeV, a 


3orn =0.30809055(l) pb 


r t = o 


u 


10" 4 


10~ 5 


10" 6 


-lloop 


0.315169(2) 


0.315167(3) 


0.315167(4) 


8 


2.297(1) 


2.297(1) 


2.297(1) 


r t 7^ 0, fixedi 


^lloop 


0.315469(2) 


0.315471(3) 


0.315472(3) 


s 


2.395(1) 


2.396(1) 


2.396(1) 



Table 4. The total lowest-order cross sections <r Born , the 
one-loop corrected cross sections <7 lloop in pb and relative 
corrections S in % at three values of Q = 10 -4 , 10 -5 , 10 -6 
and two of 





=1000 GeV, a 


3orn =0.105977185(2) pb 


n = o 


w 


10- 4 


10~ 5 


10" 6 


^lloop 


0.102872(3) 


0.102874(3) 


0.102875(4) 


S 


-2.930(2) 


-2.928(3) 


-2.927(4) 


i~t 7^ 0, fixedi 


' lloop 


0.102753(3) 


0.102775(4) 


0.102774(5) 


S 


-3.043(3) 


-3.021(3) 


-3.023(4) 



Table 5. The same but for \/I=1000 GeV. 

Tables were computed with three values of tu 
and for two options r t = 0, 7^ 0. As far as w-stability 
is concerned, we see that it depends on the /Voption. 



^=7000 GeV, a Born =2.23529503(4) fb 


r t = 




10" 5 


10~ 6 


HT 7 


^ lloop 


1.5917(1) 


1.5917(2) 


1.5919(1) 


S 


-28.79(1) 


-28.80(1) 


-28.78(1) 


-Tf 7^ 0, fixedi 


£ lloop 


1.5815(1) 


1.5895(2) 


1.5895(2) 


5 


-29.25(1) 


-28.89(1) 


-28.89(1) 



Table 6. The same but for vl=7000 GeV and for Q 
10~ 5 , 10~ 6 , 10~ 7 and the cross sections given in fb. 



For r t = there is stability for all considered y/I and 
oj. For r t 7^ it depends on the cms energy V~§: for 
200 GeV stability sets in already at Q = 10~ 4 , for 1000 
GeV at Q = 10~ 5 and for 7000 GeV at Q = 10~ 6 . 



^=200 GeV, A = 




10- 4 


10- 5 


10~ 6 




0.328077(2) 


0.328082(4) 


0.328083(4) 




6.487(1) 


6.489(1) 


6.489(1) 




0.328073(2) 


0.328078(2) 


0.328076(3) 




6.486(1) 


6.488(1) 


6.487(1) 



Table 7. The total one-loop corrected MS subtracted 
quantities a hIS in pb and corresponding <5 MS in % at Gj = 
10 -4 , 10 -5 , 10~ 6 , respectively. Subscript 1 means that light 
quark masses are used as in Eq. (|33|l . while 2 means that 
10 times smaller masses are used. 
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However, one can see that the difference between 
the two options, although it persists, is of order 1 per 
millc in absolute deviation, which is well below any 
reasonable estimate of the theoretical uncertainty. We 
therefore conclude that one may use the usual infrared 
regularization for s channel processes. 

The numbers of Tables [7]]9] were produced with the 
aim to demonstrate the independence from light quark 
masses of the subtracted quantities, a M and 6 M for 
A = option at three cms energies. One sees that the 
light quark mass independence holds in all considered 
cases at the level much lower than 1 per mille, which 
is quite sufficient for practical applications. 



^=1000 GeV, A = 




1(T 4 


io- 5 


10~ 6 


a JIs 


0.100084(3) 


0.100086(3) 


0.100087(4) 




-5.561(2) 


-5.559(3) 


-5.558(4) 




0.100084(3) 


0.100086(3) 


0.100087(4) 




-5.565(4) 


-5.567(5) 


-5.568(7) 



Table 8. The same but for V5=1000 GeV. 



4.3 t channel 

4.3.1 SANC-CompHEP comparison 

This is the most complicated case: t channel cross 
sections usually show up bad statistical convergence. 
For this comparison we use the CompHEP (v. 4. 5.1) 



Corrections to Single-top... 



^=7000 GeV, A = 




10" 5 


10~ 6 


10~ 7 




1.4566(1) 


1.4565(2) 


1.4568(1) 




-34.84(1) 


-34.84(1) 


-34.83(1) 




1.4569(2) 


1.4570(2) 


1.4570(2) 




-34.82(1) 


-34.82(1) 


-34.82(1) 



Table 9. The same but for ^=7000 GeV but for cross 
sections given in fb and for to = 10 -5 , 10 -6 , 10 -7 . 

setup, but with non-zero masses of the u and d quarks 
(accessed via be — > ts~f channel) . For the Tables of this 
subsection we used m u — = m q = 66 MeV, and 
10m q means m u = m c i = 660 MeV. The cut on the 
cms photon energy was E 1 > 2 GeV. 

Furthermore, rows marked "SANC(S)" were com- 
puted retaining m q or 10m g only in fermion propa- 
gators radiating a photon, while "SANC(F)" means 
that light quark masses were kept everywhere ( "P" ully 
massive case). Table [TT1 contains the same information 
as Table ITOl but for the process b + d — >• t + u. 

Let us discuss the results of Tables IT0UTT1 which 
are qualitatively the same, so we need not refer to the 
Table number. 

Two options (A = 0, 7^ 0) <8> (m q ) (4 upper rows) 
show very good agreement at the threshold, \/§ = 
200 GeV, and a notable deviation at TeV energies. We 
assumed that this is due to the mass singular origin 
of external light quark lines emitting photons, which 
obviously is getting worse at high energies. In this con- 
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VS/GeV 


ws 


200 


1000 


7000 


r t = 0, m q 


CompHEP 




0.2468(1) 


9.078(1) 


18.68(1) 


SANC(S) 




0.2467(1) 


9.090(1) 


18.45(1) 


r t + 0, m q 


CompHEP 


f 


0.2392(1) 


9.161(1) 


18.68(1) 


SANC(S) 


fl 


0.2392(1) 


9.174(1) 


18.70(1) 


r t = 0, 10m, 


CompHEP 




0.1740(1) 


7.512(1) 


16.11(1) 


SANC(S) 




0.1745(1) 


7.514(1) 


16.12(1) 


SANC(F) 




0.1740(1) 






r t / 0, Wm q 


CompHEP 


f 


0.1694(1) 


7.597(1) 


16.36(1) 


SANC(S) 


fl 


0.1698(1) 


7.599(1) 


16.37(1) 


SANC(F) 


fl 


0.1694(1) 







Table 10. Comparison of the cross section <7 hard (\/5, u), 
fb for the process b + u — > t + d for three cms energies; 
four options: (F t — 0, 7^ 0) <g) (m q ,10m q ); and fixed width 
scheme (ws): (f\=fixedi for the case of SMC). 



ncction, we emphasize that a special procedure of nu- 
merical stabilization of terms giving rise to light quark 
mass singularities was applied in SANC. 

To check this assumption we considered options 
{r t = 0, 7^ 0) ® (10m q ) and we immediately see the 
inverse trend: only satisfactory agreement at a/I = 
200 GeV and much better agreement at TeV energies. 
This exercise convinced us that in the mass singular 
stable regime (large masses) CompHEP and SANC results 



15 



VS/GeV 


ws 


200 


1000 


7000 


r t = 0, m q 


CompHEP 




0.1562(1) 


8.620(1) 


18.27(1) 


SANC 




0.1562(2) 


8.634(1) 


18.33(1) 


r t / 0, m q 


CompHEP 


f 


0.1522(1) 


8.706(1) 


18.53(1) 


SANC 


fl 


0.1522(1) 


8.717(1) 


18.57(1) 


r t = 0, 10m 9 


CompHEP 




0.1032(1) 


7.120(1) 


15.98(1) 


SANC 




0.1037(1) 


7.123(1) 


15.98(1) 


SANC(F) 




0.1032(1) 






r t / 0, 10m 9 


CompHEP 


f 


0.1009(1) 


7.204(1) 


16.22(1) 


SANC 


fl 


0.1014(1) 


7.207(1) 


16.24(1) 


SANC(F) 


fl 


0.1009(1) 







Table 11. The same comparison of the cross section 
<r hard (v / I, Co), fb but for the process b + d — > t + u. 



do agree. It is easy to understand why with the lat- 
ter options the situation near the threshold is not so 
ideal. In SANC we want to neglect terms with power 
dependence on light quark masses, keeping only those 
that would lead to mass singularities. Relatively large 
masses, 10m g , start to matter near the threshold. (Note 
also that on passing to the hadronic level we will need 
to subtract mass quark singularities anyway: see next 
subsection). 

To check this assumption we considered the fully 
massive variant SANC(F) with options (r t — 0, 7^ 0) <g> 
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(10m 9 ) and only at \/§ = 200 GcV. As is seen from 
the SANC(F) rows of the Tables, the good agreement 
with CompHEP is recovered. 

Although SANC(F) is identical to CompHEP, we pre- 
fer to use the much faster SANC(S) version of the code 
since the masses of u, d and s quarks are very small. 
The only exception is the c quark whose mass effect 
will be obviously notable near the threshold. Since we 
posses the SANC(F) version of the code, we can use 
a mixed variant with only c quark mass dependence 
retained. 

To conclude this section, we note that the results 
for the charge conjugate channels are identical to those 
we have considered: 





=200 GeV, a 


3orn =7.3551155(l) pb 


r t = o 


UJ 


10~ 4 


10~ 5 


10" 6 


" Hoop 
a 


7.81964(2) 


7.81954(4) 


7.81956(4) 


s 


6.3156(3) 


6.3144(5) 


6.3146(6) 


-Tf 7^ 0, fixedi 


2 Hoop 
a 


7.81954(2) 


7.81956(3) 


7.81957(4) 


s 


6.3143(3) 


6.3145(4) 


6.3147(5) 



Table 12. The total lowest-order cross sections <r Borl1 , 
the one-loop corrected cross sections <5- 11oo p in pb and 
relative one-loop correction S in % at three values of 
Q = 10 -4 , 10~ 5 , 10~ 6 and two values of A. 



b + u -> t + d to b + u -> t + d , 

b + d^t + u to b + d^t + u. (37) 

4.3.2 One-loop EW corrections, process bu — s- td 

The numerical results for this subsection were again 
produced with the SANC setup, Eq. (|33[) . with the same 
aim to demonstrate the stabilility of one-loop cor- 
rected EW cross sections a lloop and relative EWRC 
S against variation of the soft-hard separator cD and 
to study the difference between the two options: r t = 
0,^0. 

Tables [T2"Hi"4l were computed with three values of 
ui and for two options r t = 0, =/= 0. As is seen, the in- 
tervals of w-stability depend on \/§. The higher the 



^=1000 GeV, cr B 


Orn =48.99340951(6) pb 


r t = o 




10~ 5 


10~ 6 


io- 7 


^.lloop 


53.344(1). 


53.345(1). 


53.346(1). 


5 


8.880(1). 


8.881(2). 


8.883(2). 


^ ^ 0, fixedi 


^.lloop 


53.292(1) 


53.293(1) 


53.293(1) 


8 


8.773(1) 


8.776(2) 


8.776(2) 



Table 13. The same but for Vs=1000 GeV and for O 



io~ 5 ,io~ 6 ,io- 7 . 



cms energy, the lower are the values of u required 
to reach stability. For \/§=200GeV, stability starts 
at u> = 10~ 4 , for VI=1000GcV at w = 10~ 5 and 
for \/I=7000TeV it depends additionally on the r t - 
option. For r t — stability starts at Q = 10~ 6 but 
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71=7000 GeV, a Born =50.82423111(8) pb 



n = 





KT 6 


io- 7 


10~ 8 


-Hoop 


55.696(2) 


55.697(2) 


55.697(2) 


5 


9.586(3) 


9.588(4) 


9.587(4) 


r t 7^ 0, fixedi 


-Hoop 


55.613(1) 


55.639(2) 


55.640(2) 


5 


9.428(3) 


9.474(3) 


9.476(4) 



Table 14. The same but for \/5=7000 GeV and for cu = 

10~ 6 ,10" 7 ,10" 8 . 

for r t 7^ only at Oj — 10~ 7 . However, one can see 
that the difference between the two options, although 
it persists, is well below 1 per millc in absolute devi- 
ation at lower energies, reaching 1 per mille only at 
-\/5=7000 GeV, that is the difference under study is 
below any reasonable estimate of the theoretical un- 
certainty. We therefore conclude that one may use the 
usual infrared regularization also for t channel pro- 
cesses. 

Subtraction of quark mass singularities 

The numbers of the Tables [T5HT71 were produced 
with the aim to demonstrate the light quark mass in- 
dependence of the one-loop subtracted quantities a M 
and 8 M for the r t = option at three cms energies 
and three values of a). 

The stability of ct 37 ^ and S 17 ^ seen in Tables EU £7] 
is a very important check of the correctness of implc- 



^=200 GeV, a Bol ' n =7.3551155(l) pb 



r t = o 





io- 4 


10~ 5 


IO -6 




8.09370(2) 


8.09377(4) 


8.09380(4) 




10.042(1) 


10.043(1) 


10.043(1) 




8.09116(4) 


8.09127(6) 


8.09129(8) 




10.043(1) 


10.045(1) 


10.045(1) 



Table 15. The total one-loop corrected MS subtracted 
quantities (see text) a AIS in pb and corresponding 8 MS in 
% at cu = 10" 4 , 10" 5 , 10" 6 . Sub-index 1 means that light 
quark masses used are as in Eq. (|33|l . while 2 means that 
10 times lower masses are used. 



\/S=1000 GeV, cr Born =48.99340953(6) pb 


ft = 




io- 5 


10~ 6 


IO" 7 




53.416(1) 


53.416(1) 


53.417(1) 


5 MS 


9.027(1) 


9.028(2) 


9.030(2) 




53.416(1) 


53.416(1) 


53.417(2) 




9.029(2) 


9.029(2) 


9.032(3) 



Table 16. The same but for \/I=1000 GeV and for Q = 

io _5 ,io- 6 ,io- 7 . 

mentation of the complete EWRC and of the initial 
quark mass singularity subtraction. The latter is cru- 
cial at the envisaged stage of going to the hadron level, 
i.e. convolution with parton density functions. 
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71=7000 GeV, a Born =50.82423111(8) pb 



n = 





10" 6 


io- 7 


10 -8 


& MS 


55.699(2) 


55.700(2) 


55.699(2) 


5p 


9.591(3) 


9.593(4) 


9.592(4) 


a?* 


55.698(2) 


55.697(3) 


55.696(3) 


sW 


9.589(3) 


9.588(5) 


9.585(6) 



Table 17. The same but for 71= 7000 GeV and for Q = 

10 _ 6 ,10 76 ,10 -8 . 

4.3.3 One-loop EW corrections: process bd —> tu 

Here we repeat the study of subsection 14.3.21 but now 
for the crossed t channel. Tables [T5H2H1 are analogu- 
ous to Tables [T2T[T4"1 They demonstrate the stabilility 
of a 11 and o~ against variation of U) and show the 
difference between the options r t = 0, 7^ 0. 

Tables ITSl - f^Ul show up properties very similar to 
the Tables for the direct t channel. Each Table shows 
the presence of a plateau of stability in u and the 
same level of difference between calculations for two 
options r t = 0,7^ 0, qualitatively with very similar 
energy dependence. 

Subtraction of quark mass singularities 

The numbers of Tables I2TH231 were also produced 
with the aim to demonstrate the light quark mass in- 
dependence of the subtracted quantities a M and i5 M 
for option r t = at three cms energies. 



71=200 GeV, cr Born =4.495790646(3) pb 



r t = 





10- 4 


10~ 5 


10" 6 


2 Hoop 
a 


4.87169(2) 


4.87165(2) 


4.87165(3) 


S 


8.3612(4) 


8.3603(5) 


8.3603(7) 


i~f 7^ 0, fixedi 


2 Hoop 
a 


4.86989(5) 


4.86990(2) 


4.86991(3) 


s 


8.3212(4) 


8.3213(5) 


8.3215(6) 



Table 18. The total lowest-order cross sections <r Born , 
the one-loop corrected cross sections <5- 11oo p in pb and 
relative one-loop correction S in % at three values of 
w = 10 -4 , 10~ 5 , 10~ 6 and two values of A. 



7§=1000 GeV, cr B 


3rn =46.68695597(6) pb 


r t = 




10~ 5 


10 -6 


IO" 7 


^.lloop 


51.054(1) 


51.054(1) 


51.055(1) 


6 


9.354(1) 


9.354(2) 


9.356(2) 


i~f 7^ 0, fixedi 


^.lloop 


51.004(1) 


51.005(1) 


51.006(1) 


8 


9.246(1) 


9.249(2) 


9.250(2) 



Table 19. The same but for 71=1000 GeV and for lo = 

10- 5 ,10" 6 ,10" 7 . 

Again, Tables l2"TH2"3l although within larger statis- 
tical errors, show the stability of a M and 5 M for all 
cms energies and Q values considered, confirming the 
correctness of implementation of the complete EWRC 
also for the case of the crossed t channel. 
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71=7000 GeV, a Born =50.72449055(7) pb 



n = 





KT 6 


io- 7 


10~ 8 


~ lloop 


55.581(2) 


55.583(2) 


55.580(2) 


5 


9.575(3) 


9.578(4) 


9.572(4) 


ft 7^ 0, fixedi 


-Hoop 


55.495(2) 


55.527(2) 


55.528(2) 


5 


9.404(3) 


9.467(3) 


9.469(4) 



Table 20. The same but for 75=7000 GeV and for cu = 

10~ 6 ,10" 7 ,10" 8 . 



71=200 GeV, a Born =4.495790646(3) pb 


n = 




i<r 4 


nr 5 


IO -6 




4.92749(2) 


4.92749(2) 


4.92749(3) 


S MS 


9.6023(4) 


9.6023(5) 


9.6023(7) 




4.92515(3) 


4.92513(4) 


4.92514(5) 




9.6065(6) 


9.6061(9) 


9.606(1) 



Table 21. The total one-loop corrected MS subtracted 
quantities a MS in pb and corresponding S MS in % at 
Q — 10 -4 , 10 -5 , 10 -6 . Sub- index 1 means that light quark 
masses used are as in Eq. (|33|) . while 2 means that 10 times 
lower masses are used. 

5 Conclusions and Outlook 

In this paper we describe the implementation into the 
SANC framework of the complete one-loop EW calcu- 
lations, including hard bremsstrahlung contributions, 
for the processes of the top quark decays a of s and 



71=1000 GeV, cr Born =46.68695597(6) pb 



r t = 





10~ 5 


10~ 6 


io- 7 




51.086(6) 


51.0858(8) 


51.087(1) 


sW 


9.422(1) 


9.422(2) 


9.424(2) 




51.0856(9) 


51.086(1) 


51.087(2) 


sW 


9.424(2) 


9.424(3) 


9.428(3) 


22. The same but for 75=1000 GeV and f 
0- 6 ,hT 7 . 


75=7000 GeV, a Born =50.72449055(7) pb 


ft = 




nr 6 


io- 7 


10 -8 


-MS 


55.583(2) 


55.585(2) 


55.582(2) 


5 MS 


9.579(3) 


9.582(4) 


9.575(4) 


&F 


55.583(2) 


55.585(3) 


55.576(3) 


5p 


9.579(4) 


9.582(5) 


9.564(7) 



Table 23. The same but for 75=7000 GeV and for Q = 
10~ 6 , 10 -7 , 10~ 8 . 

t channel production, the latter two at the partonic 
level,. The essentially new aspect of this paper is the 
study of regularisation of the top quark-photon in- 
frared divergences with aid of the complex mass of 
the top quark. For this reason, we briefly come back 
to the top decay channels, since this issue was not 
studied in our previous papers, Rcf. [TJH]. A compari- 
son of these electroweak corrections computed within 
this new approach with those computed by the con- 
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ventional method showed a sizable (1%) effect for top calculations were done in gauge), checked the sta- 

decays and only per mille level effect for both consid- bility of the result against variation of the soft-hard 

ered top production processes. In this paper we limit separation parameter u and the independence of the 

ourselves to 4-fermion type of processes. For this rea- MS subtracted quantities off the initial quark masses 

son we do not consider t quark production channel in which is crucial for calculations at the hadronic level, 
association with W boson. This channel will be pre- 



As has become SANC standard, we tried to compare 
our numerical results for these channels with other 
independent calculations. For the decay channels it 
was done in our preious papers (both EW [Tj and 
QCD [53]), showing good agreement. As usual, the 
Born level and the hard photon contributions of all 
three channels were checked against CompHEP package 
and we found very good agreement for both _T f = 0, 7^ 
options. 



scntcd elsewhere. 

We have presented analytical expressions for the 
covariant amplitude of the process and the helicity am- 
plitudes for three different cross channels within the 
standard SANC multichannel approach, Ref [17] . In this 
approach the one-loop covariant amplitude is com- 
puted only once for to-the-vacuum-annihilation pro- 
cesses, see Figs. [T] and To get the CA for a physical 
channel a cross transformation of external 4-momenta 
is performed. The helicity amplitudes have been cal- 
culated for each channel separately. 

Within the SANC framework, we have created the 
standard FORM and FORTRAN modules, see, Ref. 
|22| . compiled into a package sanc_cc_vl .40 which 
may be downloaded from SANC project homepages. All 
the calculations were done using a combination of ana- 
lytic and Monte Carlo integration methods which will 
make it easy to impose experimental cuts in furtcom- 
ing calculations for the pp-collisions at LHC, which 
will be the subject of a subsequent paper. 

The emphasis of this paper is to be assured of the 
correctness of our results. We observe the indepen- The results presented in this papers lay a solid base 



As far as the comparison of EWRC for the pro- 
duction processes is concerned, the situation is not 
satisfactory. We did not find in the literature any cal- 
culations for the s channel. For the t channel there 
are many papers, see, for instance, Refs. [121(15] and 
references therein. The most appropriate is Ref. |12j . 
where however the hard photon contribution is not 
included. We tried to suppress it in our calculations, 
but got only rough agreement with their Fig. 3(a). The 
most advanced paper Ref. [Tl] contains results at the 
hadronic level and, at present, we can not yet make a 
comparison. 



dence of the form factors on gauge parameters (all for subsequent extensions of calculations for the single 
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top production channels at hadron colliders, basically 10. S. Frixione, E. Laenen, P. Motylinski, and B. R. Web- 
LHC. ber, JHEP 03 (2006) 092, |hep-ph/0512250| 
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